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ABSTRACT 

There is still no consensus as to what causes galactic discs to become warped. Suc- 
cessful models should account for the frequent occurrence of warps in quite isolated 
galaxies, their amplitude as well as the observed azimuthal and vertical distribution 
of the Hi layer. Intergalactic accretion flows and intergalactic magnetic fields may 
bend the outer parts of spiral galaxies. In this paper we consider the viability of these 
non-gravitational torques to take the gas off the plane. We show that magnetically gen- 
erated warps are clearly flawed because they would wrap up into a spiral in less than 
two or three galactic rotations. The inclusion of any magnetic diffusivity to dilute the 
wrapping effect, causes the amplitude of the warp to damp. We also consider observa- 
tional consequences of the accretion of an intergalactic plane-parallel flow at infinity. 
We have computed the amplitude and warp asymmetry in the accretion model, for a 
disc embedded in a flattened dark matter halo, including self-consistently the contri- 
bution of the modes with azimuthal wavenumbers m = and m = 1. Since the m — 
component, giving the U-shaped profile, is not negligible compared to the m = 1 com- 
ponent, this model predicts quite asymmetric warps, maximum gas displacements on 
the two sides in the ratio 3 : 2 for the preferred Galactic parameters, and the pres- 
ence of a fraction ~ 3.5% of U-shaped warps, at least. The azimuthal dependence of 
the moment transfer by the ram pressure would produce a strong asymmetry in the 
thickness of the Hi layer and asymmetric density distributions in z, in conflict with 
observational data for the warp in our Galaxy and in external galaxies. The amount 
of accretion that is required to explain the Galactic warp would give gas scaleheights 
in the far outer disc that are too small. We conclude that accretion of a flow with no 
net angular momentum, cannot be the main and only cause of warps. 

Key words: galaxies: intergalactic medium - galaxies: kinematics and dynamics - 
galaxies: spiral - galaxies: structure - magnetic fields 



1 INTRODUCTION 

The generation and persistence of galactic warps have been a puzzle for galactic dynamicists over the last decades (see 
Toomre 1983; Binney 1991, 1992 for reviews) and it is still a problem of huge activity. Intergalactic winds, tidal interactions 
with satellites, magnetic fields, misaligned dark matter haloes, cosmic infall, vertical resonances and bending instabilities 
have all been proposed as causative agents (e.g., Binney 1992; Griv et al. 2002; Revaz & Pfenniger 2004). Understanding 
the underlying physics responsible for the generation of galactic warps may shed light on the late stages of accretion in 
galaxies or the interaction between the periphery of galaxies and the teneous intergalactic medium. Motivated by some recent 
observational correlations that favour scenarios in which the intergalactic medium is a necessary ingredient to understand 
the warp mechanism (Castro- Rodriguez et al. 2002; Sanchez-Saavedra et al. 2003), we have reconsidered the role of the 
intergalactic medium. 
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Figure 1. Schematic diagram showing the shape of the warp for the wind model and for the accretion model. 

The idea that the Galactic gas layer bends by a subsonic (intergalactic) wind was proposed by Kahn & Woltjer (1959). 
In this picture, which we refer to as the wind model, the discs of the galaxies are distorted because of pressure differences 
exerted by a wind flowing past the "solid" halo of hot gas. No accretion of mass in the disc occurs in this model because the 
coronal gas acts as an obstacle. Saar (1979) shown that if the coronal gas responses against the colliding wind as a smooth 
halo, which is a better representation than a rigid obstacle, the warping effect is considerably reduced. A major problem with 
the wind model is that the required massive wind, for which there is no evidence, should be observable (Binney 2000). 

Lopez-Corredoira et al. (2002), hereafter referred to as LC, have explored in some detail the idea proposed by Mayor & 
Vigroux (1981) that the accretion of baryonic matter, probably in the form of High Velocity Clouds (HVCs), as the galaxy 
moves through the intergalactic medium, can explain the existence of warps. Unlike the wind model, the intergalactic gas 
flow is supposed to be accreted directly onto the disc creating a "coUisional torque" . We will refer to the latter scenario as 
the accretion model. 

The wind and accretion models have in common that the impulse transfer in an oblique incident flow, has an azimuthal 
dependence and, hence, it can act as a source of excitation of warps. However, the azimuthal dependence of the impulse 
transfer is not the same for both models. For a flow oriented as in Fig. 1, the wind model predicts that the left side of 
the galaxy would lie above the plane defined by the inner galaxy and the right side should be placed below this plane (see, 
e.g., Ikeuchi & Tomisaka 1981). However, in the accretion model, it holds the opposite. Therefore, since their effects should 
partially counteract, the amplitude of the warp depends on what fraction of the intergalactic gas resides in clouds that can 
reach the disc without evaporating. 

It is also important to remind that the accretion model is absolutely different to the cosmic infall one (e.g. Ostriker & 
Binney 1989). These authors pointed out that the accretion of spherical shells of protogalactic material in the outer halo 
shifts significantly the direction of the net angular-momentum vector of a typical galactic halo every Hubble time and that 
the resulting gravitational torques are liable to warp the disc. 

An alternative model that invokes the intergalactic medium as a potential source to induce warps was provided by 
Battaner et al. (1990). These authors suggested that the tension of galactic magnetic fields that arises in the periphery of 
galactic gas discs when the field lines exit the disc and connect with the uniform intergalactic field, produces the required 
torque. This magnetic model predicts that the gas layer should bend more strongly than the stellar disc, which motivated 
Cox et al. (1996) to observe the galaxy UGC 7170 in both optical and radio (see also Casertano et al. 1987). They did not 
see any noticeable difference between them, but a single object is not enough to rule out the magnetic hypothesis in general. 
On the other hand, Florido et al. (1991) found a clear wavelength dependence in the warp's curve for three galaxies. This fact 
could suggest that warps are induced by a force acting on the gas but not on stars. 

A virtue attributed to these models (the magnetic scenario, the wind and the accretion model) is that they can explain 
the alignment of waxps of different galaxies in a natural way (e.g., Battaner et al. 1991). Other advantage of these models is 
that they work with no need for the existence of a massive dark halo of coUisionless particles (see LC; Castro-Rodrfguez et 
al. 2002). 

The first question that arises is whether waxps trace fairly the cosmic infall or, on the contrary, intergalactic accretion 
flows and magnetic forces arc working sirnultanoously, either coherently or destructively, as bending mechanisms. In the latter 
case it would be interesting to constrain the contribution of each agent. It seems unlikely that all the different warping 
mechanisms may be operating simultaneously with comparable strength. It is important to ascertain the role of the different 
torques in modelling the outer parts of disc galaxies, being aware of the difficulties inherent to each scenario. In this paper we 
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will analyse the eflFects of the non-gravitational torques. In §2 we provide the departure point for our study of the steady-steate 
configurations of gas discs under the action of torques. Section 3 is devoted to the magnetic model, whereas the accretion 
model is dcfored to §4. Wc conclude that the effect of non-gravitational torques, either magnetic or hydrodynamical by the 
intergalactic medium, is not the main cause for warps generation in field large spirals. Exceptionally, cluster galaxies suffering 
from strong ram pressure stripping could be bent by dynamical pressure. 



2 TORQUES ON GALACTIC DISCS. BASIC CONSIDERATIONS 

Kahn & Woltjer (1959) shown that, as galactic warps cannot be primordial because of the so-called winding problem, an 

external torque to take the gas off the plane should be required. These torques can be due to tidal forces, magnetic fields, by 
ram pressure or by cosmic infall of matter in the outer halo. We must warn that the presence of torques is not guarantee of 
the existence of an immortal steady antisymmetric warp. It is essential to study the evolution of the disc to ensure that the 
initial forced integral-sign form of the warp is not dispersed. In this Section we will discuss some basic points on the evolution 
of a gas disc subject to torques, which are necessary to understand the remainder of the paper. 



2.1 External torques and coordinate system 

The variation of the angular momentum of the whole galactic disc idisc (stars plus gas) will depend only on the external 
torques according to the equation: 

= Th + Tp + r„g, (1) 

where th is the integrated torque over the disc produced by the dark matter halo, Tp the gravitational torque due to the 
presence of a perturber (as a companion galaxy) and Tng is the torque caused by a non-gravitational force, /^g, which may 
have a magnetic origin (see §3) or may describe the action of the dynamical pressure applied by the intergalactic medium onto 
the disc (§4). In the present work we investigate the role of the force f^^ in shaping the warp of galaxies in tidal isolation 
such that /p = Tp = is a good approximation. 

For a galaxy embedded in a spherical halo we have th = 0, hence 

dt ~ 

So, in a stationary configuration, the disc must precess with a no-null frequency u). Nevertheless, for the models considered 
in this paper, u) is very small (see §4.1 and LC). In particular, the frequency u! in both the accretion and magnetic models is 
of 1 cycle per 32 Gyr. This period is much larger than the characteristic dynamical time-scale in the disc, < 1 Gyr, and 
refiects the fact the the disc is very reluctant to precess as a whole. Once the precession of the disc is neglected, the inner 
or central disc (up to ~ 10-12 kpc in the Galaxy), which is almost unaffected by the external torque, defines our inertial 
coordinate system, being the z-axis perpendicular to the central disc. The error made by this approximation in the calculation 
of torques in the warp is 0(w/n), where Q is the angular speed. This approximation allows us to take the inner disc as being 
non-responsive so that the gravitational force due to the inner disc is incorporated in the model as an external background 
force acting on the outer disc, and consider only the evolution of the the outer parts of the galaxy, where the surface density 
of the disc drops quickly with radius and the disc starts to bend. The equation for the angular momentum of the outer disc, 
is 



dL 



(out) 

disc (out) , (out) 



dt 



■ T 



+ t\Z\ (3) 



where Tj°^*^ is the torque produced by the inner disc onto the outer disc. If, for some reason, it holds rig"*^ = — tJ°|J'', then 
the mean precession of the outer disc will be zero. However, if we model the outer disc as a collection of concentric rings and 
since galactic discs do not behave as rigid solids, each ring may precess at a different frequency, even though the external 
torques cancel out. 

On the other hand, if the dark halo is oblate or prolate and if it is taken to be frozen so that the effect of the disc gravity 
on the halo is thus neglected, a situation th = —Tng is possible, implying that dLdiec/dt = 0. In that case the mean precession 
of the outer disc vanishes, as external torques are in balance. 



2.2 Equation of motion and bending modes 

It is convenient to make discussions in terms of the vertical displacement of the disc at {R,4},t), denoted by h{R,4i,t), where 
R is the perpendicular distance from the a-axis and (f) is the azimuthal angle measured in the equatorial plane. As it is usual 
(e.g. Hunter & Toomre 1969; Sparke & Casertano 1984), we will assume throughout this paper that the vertical scale is 
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large enough compared to the horizontal epicycle that the epicyclic motion in R and 4> can be neglected. The entire vertical 
dynamics can be summed up by the equation 

^ + n^) h = Ui + U + U^, (4) 

where fseii{R,4'jt) and fh{R,(l>,t) represent gravitational vertical forces per unit mass due to the self-interaction of the disc 
material and interaction with the dark halo, respectively. fng{R, (f>,t) is a non-gravitational force (per unit mass) to be specified, 
and Q{R) is the angular speed about the axis i? = 0. In practice, since the disc is assumed to stay close to the midplane z = 0, 
i.e. h <^ R, /ng does not depend on h to first order. 

To first order in h/R, the vertical force f sell {R,(t>,t) due to the disc itself is given by 

where A{R, R' , 0, 0') = {R^ + R''^ - 2RR' cos((^ - </)') + zlf^ and zo is a softening radius that arises because the disc has 
finite vertical thickness. If the dark halo is axially symmetric about R = Q and mirror symmetric about the plane z = 0, fh 
can be approximated by 

fh{RA,t) = -vlh{RA,t), (6) 

where Vh{R) = d^<^h/dz^ evaluated at the equatorial plane, with ^h{R, z) the gravitational potential created by the dark 
halo. After specifying /ng and the initial conditions of the disc, the subsequent motion of the disc is completely described by 
Eq. (4). 

First, we will focus on modes with azimuthal wavenumber m = 1 of the form h{R, (f>, t) — hi{R) sin(0 — uiRt + rj), where 
the eigenfrequency lvr may depend on R, and ?7 is a phase which could depend also on R but it is chosen constant in order to 
have an integral sign warp at t = 0. Each mode consists of a collection of concentric tilt rings that precess with frequency u)r. 
If the frequencies of precession uir vary between successive rings, the line of nodes, where material at each radius encounters 
the reference galactic plane, will wrap in a spiral form in a characteristic time ~ n/[u)R{Ri) — ujr{R2)], where Ri and R2 are 
the minimum and maximum radii of the warp. In order to avoid the winding problem, we are mainly interested in modes in 
which the galaxy can process as a reasonably rigid unit or, conversely, that ujr is approximately constant in R. 

If /ng = 0, the gravitational torques exerted in each ring drive them to precess with uir ^ about ii = 0. It is 
straightforward to estimate in the simplest approximation that the gas in the outer disc moves as test particles in a fixed 
axisymmetic galactic potential, then /self = —Vdh- Equation (4) is simplified to 



where vl{R) = v^{R) + vl^{R). After direct substitution of h into Eq. (7), we find that the rings precess with eigenfrequencies 
ijjR = f2 ± z^t. For an oblate potential vt > O,, so the rings regress with angular frequency i>t — Q, or advance with frequency 
Vt + Kahn & Woltjcr (1959) soon noticed that for a galaxy like the Milky Way, the warp will wind up in a characteristic 
time < 2 Gyr, which is very short compared to the age of the disc (see also Binney 1992). 

Sparke & Casertano (1988) shown that if one uses Eq. (5) instead of the approximation fseit = — i^Jft, a self-gravitating 
disc embedded in a fiattcncd dark halo can support a discrete tilted mode since it is able to reconcile the divergent precession 
rates of inner and outer rings. However, Nelson & Tremaine (1995) found that when the response of the dark halo to the 
precession of the disc is included, the modified tilted mode is damped within a few dynamical times (see also Binney et 
al. 1998). 

In the next Subsection, we will study the bending waves of a disc which is subject to an extra force /ng. The existence 
of modes in which the disc can precess as a whole depends on the form of /ng. As /ng does not necessarily depend linearly 
on h, wo refer to /ng as the "forcing term" of the complete inhomogcncous equation (4). When /ng 7^ 0, the solution will be 
a superposition of the above-mentioned modes arising from the homogenous equation (/ng = 0) plus the particular solution. 
Since those waves associated with the homogeneous equation are unable to survive as long-lived configurations, we will be 
only interested in the particular solution. 



2.3 Stationary conflgurations 

For reasons that will become clear later, we will consider two cases: case A and case B. In case A, we assume that the force 
/ng does not depend on time, i.e. /ng = fiig{R,(f>). Since we wish to study only distortions with wavenumber m = 1, it is 
natural to restrict the present analysis to a vertical force with m — 1. Without the loss of generality, this is accomplished by 
a force of the form /ng(-R, 4>) = fi{R) sincj), where fi{R) represents an arbitrary function of R. Case B assumes that the force 
is described by fng{R,4>,t) = /i(i?)sin(0 — Q,t). Case A and B are relevant for the wind and magnetic models, respectively 
(see §§3 and 4). 
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Case A.- Adding the vertical force, the dynamical equation for h is 

/l = -^/t'/^ + /l(i^)sin<A. (8) 

Here, for the sake of clarity, we have assumed that /sdf = —v^h (see §4 for a more realistic approximation for /self)- The 
solution of this equation will be the sum of the solution of the homogeneous equation (Eq. 7) , plus the particular solution hp 
given by 

K{RA,t) = ^^s,m^, {ut^n). (9) 

Restricting ourselves to modes with azimuthal wavenumber m = 1, there are three azimuthal frequencies for vertical excita- 
tions: vt — ^i, i^t + n and (always present). For arbitrary initial conditions the evolution is more complex. For illustration, 
in Appendix A we provide the evolution of a disc that initially lies in the equatorial plane with zero vertical velocity. As 
already said, we will retain only the particular solution as those bending waves associated with the homogenous equation will 
disappear rapidly as a result of the winding process aided by dissipation due to cloud-cloud collisions ^ . 

The particular solution, which is independent of the initial conditions and the only observable distortion at later times, 
is a zero-frequency mode (rings do not precess at all) with a straight line of nodes aligned with the x axis. The orbits tilt 
around the y axis towards the x axis only close to the resonance vt ~ condition which is not met in galactic discs. For 
a description of the particular solution in terms of orbits, we refer the reader to Appendix A. Essentially, a cloud under the 
action of a periodic vertical force executes a motion along z which is a combination of two oscillations, one with the intrinsic 
frequency ut and one with frequency O of the force /ng. All wc say is that if the amplitude of the intrinsic mode is zero, the 
cloud moves with the same frequency Q in the three directions x, y and z, thereby describing a closed circular orbit in a plane 
inclined to its galaxy's equatorial plane and rotated around the x-axis. 

The stationary solution h{R,(f>) = {i^i — ^^)^^ fi{R) sintf) describes a collection of tilt rings which do not precess. One 
could wrongly think that a situation with ujr = is not possible because each ring is subject to an external torque and hence 
all the rings must necessarily precess. There is no any inconsistency; the total (disc+halo+ applied by fng) torque acting on 
each ring is zero in this particular configuration (see Appendix A for details). In fact, integrated over a ring, the force /ng 
produces a torque directed around the a;-axis, but this torque is balanced by the torque created by the gravitational restoring 
forces /self + fh- 

In conclusion, for a forcing term of the type /ng {R, 4>) — fi (R) sin the only long-lived warped configuration has a zero 
precession rate. 

Case B.- We consider now a case in which the non-gravitational force can be expressed as f{R,4>,t) = fi{R) sm{4>— Qt). 
With the approximation /self = —Vdh, the corresponding equation for h is 

§1+^-^^ h = -ly^ h + MR) sin {4>-nt). (10) 

The non-gravitational force that feels a certain cloud is independent of time since both the cloud and the force rotate around 
the «-axis with the same angular frequency f2. 

Following the same reasoning than in case A, only the particular solution given by 

hp{R,4>,t) = sin {(I,- nt), (11) 

is considered. The corresponding vertical velocity of the fluid for this solution is 

^.+^4l]h = 0- (12) 



^dt d(f>j 

According to Eq. (11), the vectors normal to the rings precess about the z axis with frequency r2(_R) and, since the frequencies 
of these modes differ for successive rings in a disc with differential rotation, their relative phases will evolve in time and the 
bending packet will disperse in a characteristic time ~ 7r/[n(i?i) — Q.{R2)], where R\ and R2 are the minimum and maximum 
radii of the warp. For a typical galactic disc with a flat rotation curve, this time is oven shorter than that for those modes 
arising from the homogeneous solution. Consequently, this form of forcing cannot handle a long-lived warp. 

Sparke & Casertano (1988) shown that an isolated (i.e. /ng = 0) self-gravitating disc embedded in a flattened halo can 
support a discrete tilted mode since it is able to reconcile the divergent precession rates of inner and outer rings. If hsc {R, 4'^ ^) 
is the Sparke & Casertano normal mode of oscillation with /ng = 0, then the solution of the dynamical equation for m = 1 
waves, after adding the forcing term, has the form 



^ In other words, we are not interested in the initial-value problem but in the existence of a mode in which the galaxy can precess as a 
reasonably rigid unit. This is very common in other areas of physics as when deaUng with the propagator of particles in the Feynman 
diagrams, where a pole is introduced in order to have a well-defined momentum before and after interactions. 
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h{R, 4>, t) = /isc (i?, (t>, t) + h'p{R) sm{<j) - Qt), 
where hp satisfies the integral equation 

fi{R) - vlh'p{R) = Gh'p{R) / J:{R')H{R,R')R'dR' 



(13) 







G / S{R')J{R,R')h'p{R')R'dR', 



(14) 



with 



H{R,R') = [ 



(i?2 + - 2RR' cos (j) + 0g)3/2 ' 



d(t> 



(15) 



J{R,R')= f 



(ii2 + - 2RR' cos c;* + ag)3/2 " 



cos (f> d(f) 



(16) 



The discrete normal mode hsc will be damped within one dynamical time of the disc by dynamical friction with halo particles 
(Nelson & Tremaine 1995), whereas the bending waves due to forcing, being of the form sm{(j) — Q.t) with 0(7i) oc 1/i?, will 
also dissapear quickly by the wrapping process. Therefore, even including self-gravity, a torque of type B cannot contribute 
to alleviate the problem of the persistence of warps. 



3 GALACTIC WARPS IN THE MAGNETIC SCENARIO 

Battaner et al. (1990) pointed out that extragalactic field lines entering the outer parts of the gas discs could lead to a tension 
able to produce a torque. The attractive points of this model have been stressed by different authors (Battaner et al. 1991; 
Kuijken 1991; LC; Castro-Rodn'guez et al. 2002). If this mechanism is at work, waxps could trace out the direction of ordered 
extragalactic magnetic fields at Mpc scales (e.g., Vallee 1991). LC argue that a field strength of ~ 1.4 jxG could explain 
the warp's curve of the Milky Way very well and, therefore, this mechanism should not be taken lightly. Since observational 
determinations of the strength of intergalactic magnetic fields are very difficult and scarce, it would be desirable to make more 
theoretical predictions and compare with observations to quantify how much magnetic fields can contribute to support warps. 
In Section 3.1, we will show that magnetically driven m = 1 distortions will wrap in a spiral in a short time-scale due to the 
stretching by differential rotation. Turbulent or ambipolar diffusion can be invoked to dilute the undesirable wrapping effect. 
In that case, however, the amplitude of the warp will decay in a shorter time-scale (§3.2). Before doing so, we will briefly 
review the hypothesis behind the magnetic scenario for warps formation. 

For our purposes we only need to consider the vertical component of the equation of motion of a gas disc embedded in a 
general large-scale magnetic fleld: 



Here d/dt is the covariant derivative and P = Pt -\- Pm, with Pt the gas pressure consisting of the turbulent and thermal 
components and Pm is the magnetic pressure. Consider now the third term of the RHS of Eq. (17). The total vertical force at 
{R, 4>), due to this term, integrated along the vertical direction of the layer between z = —zo and zo, with zo much larger than 
the typical scaleheight of the disc, should be very small or zero as Jll^ dB^ jdz dz ~ for a disc embedded in a homogeneously 
magnetized ambient medium. Therefore, this term could contribute to flare the disc but is unable to generate either a net 
force or a torque in the disc. A similar argument holds for the pressure gradient term. Hence, integrated over all z, the net 
vertical acceleration of the disc is 



where p = with S(7?) the disc's surface density of gas and H{R) the thickness of the disc. For simplicity, we have 

adopted the simplest assumption that the potential created by the disc and the dark halo is fixed. 

In the magnetic model for the formation of warps, the gas disc is bent in an integral-sign warp by the term BbOBz / dR 
(Battaner et al. 1990). The reason is as follows. The extragalactic magnetic field in the open space, denoted by Be, is assumed 
to be uniform and plane-paxallel, forming an angle a with the plane of the galactic central disc. Without loss of generality, 
suppose for convenience that B^.z > 0. Within the thin, inner disc (say 7? < -Ri), the field is taken as approximately azimuthal, 
with Bz « 0, as observations suggest, and, therefore, the magnetic force terms in Eq. (18) equal zero. In the outer parts of 
the disc, the galactic field lines should connect with the ambient intergalactic magnetic field. In order to match the field 
configuration in the inner disc where w into the extragalactic magnetic fleld Be,z, which is nonvanishing whenever a is 




(17) 




(18) 
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(a) (b) 

Figure 2. Vertical collapse of a protogalaxy as an illustration of the formation of the warp in the magnetic scenario, (a) The protogalaxy 

before collapse, showing a uniform intergalactic field threading it. The magnetic field makes an angle a with the rotation axis Q. (b) The 
magnetic field lines (solid line), together with the warp, after collapse to the galactic disc. A more detailed description of intermediate 
stages can be found in Ifoward & Kulsrud (1997). During the time when the galaxy contracts we ignore any rotation. The turbulent 
diffusion acts to straighten out the magnetic field lines (dashed line), leading to the reduction of the magnetic tension (see §3.2). 



not close to or to 7r/2, the magnitude of the z-component of the magnetic field should increase radially outwards until it 

connects with the extragalactic magnetic field. Therefore dBz/dR is positive within Ri < R < R2, where R2 is the radius 
of the end of the disc. As is a consequence of the penetration of the extragalactic field, which is supposed to depend 
only on R, is taken independent of 0. For the same reason Bji{R, 0) = Bjifi{R) + Bji^\{R) sin[^ + with C, an arbitrary 
constant phase. Since dB^/dR > 0, the term BudB^/dR produces a torque in the outer parts of the discs, whereas the term 
B^dBz/dtf) is supposed to be negligible. Since the m = mode is not interesting as far as integral-sign shape warps concern, 
we will adopt the simplest assumption Br^ = 0. In summary, the initial magnetic field (at t = 0) in the tilted disc takes the 
form: 

B{R, (j,,t = 0)= Br,i{R) sin(0 + C)eR + B^R, ^)e^ + B4R)e^, (19) 

with dBz/dR > 0, and B when R 00. This initial magnetic configuration naturally arises in the vertical collapse 

of a spherical protogalactic cloud with a uniform magnetic field threading it (see Fig. 2), if one ignores any rotation. A more 
detailed description of this process can be found in Howard & Kulsrud (1997) and there is no need to be repeated here. The 
resulting magnetic field together with the warp are sketched in Fig. 2b. Therefore, in order to have the intergalactic field 
penetrating into the H I disc, we do not have to invoke any magnetic diffusion to transport magnetic field from outside space 
into the disc. 

As we have anticipated in §2, the presence of torques is not a sufficient condition for the warp to rotate as a rigid unit. In 
order to assess the viability of this model, we need to study the temporal evolution including the gravitational and magnetic 
torques. In the original treatment of Battaner and collaborators (also in LC), they assume that the magnetic configuration 
is unaffected by differential rotation of the disc and ignore any magnetic diffusivity as well as the vertical motions associated 
with the warp, so that dB/dt = 0. In the following, we will attempt to ascertain the consequences of the inclusion of the 
motions of the gas. These motions alter the magnetic field lines, and therefore, the form given in Eq. (19) holds only at t — 0. 

The general problem of the evolution of a magnetized rotating disc is extremely complex. Nevertheless, as illustrated 
in §2, when analyzing the success or failure of an external torque to conduct warps, it is assumed that the solution of the 
homogeneous equation gets wrapped so rapidly that in practice only the particular solution could be observed. The case of a 
magnetic torque differs from the examples in §2 in the fact that the magnetic field and the velocity field are not independent 
variables but linked through the induction equation. In other words, the force depends also on the velocity field of the fluid. 
In Appendix B we show that, in the magnetic case, vertical motions can distort the magnetic configuration, spoiling out the 
initial sirmsoidal-form (m = 1) of the vertical magnetic force. The aim of this Section is not to study the evolution of a disc 
starting with a warp of arbitrary amplitude but to find a particular solution. In fact, if we start with a force with azimuthal 
wavcnumber m = 1. we seek the particular solution that docs not generate harmonics of higher wavcmimber. In Appendix B 
it is found that the particular solution corresponds to a case where the vertical velocities are very small or zero. Thus, we will 
adopt Vz = in the next Subsections. The characteristic evolutionary time-scales for the magnetic field when Vz =^ are also 
given in Appendix B. 



8 F. J. Sanchez- Salcedo 



3.1 Differential rotation 

We first focus on tlie effect of pure differential rotation and ignore for a moment otiier effects sudi as turbulent diffusion. As 
it will be discussed in §3.2, if turbulent diffusivity is unquenched, the role of turbulent motions may also have a significant 
influence on the evolution of the warp's amplitude. Suppose that a disc was initially forced to warp by a magnetic torque so 
that it presents a straight line of nodes. We wish to know whether the line of nodes keeps straight or, on the contrary, gets 
a spiral form as time goes by. For this purpose we need to include the evolution of the magnetic field. For frozen-in magnetic 
field lines, the magnetic field evolves according to 
FiTi 

^ = Vx{vxB), (20) 

where is the ordered velocity of the gas, v — v^e^ + v^e^, and any contribution to the large-scale magnetic field by turbulent 
motions has been ignored as a first approximation. According to the ongoing discussion, we seek the particular solution with 
Vz = 0. The equation of magnetic induction with v = V4,e^ = RD.e^, in components is 

To derive Eqs. (21)-(23), it was also assumed for simplicity that the galaxies follow cylindrical rotation, i.e. dv^jdz ~ 0. 
Following a fiuid element that moves with the galactic rotation, only the azimuthal component B^j, will increase due to the 
stretching of the radial magnetic field by the differential rotation until it reaches saturation by the combined effects of 
ambipolar and turbulent diffusion. After integration of Eqs. (21) and (23) with the initial conditions Br = -Bi{,i(i?) sin(0 + C) 
and Bz{R) a function only of R, and substituting into Eq. (18) with dBz/dcj) = 0, the equation of motion reads: 

| + n|)\ = -.,^.+ |M^,in(,_a. + c). (24) 

By comparing Eq. (24) with Eq. (10), one immediately realizes that the magnetic torque corresponds to case B. The particular 
solution obeys 

hp{R,4>,t) = ^^.^sm{<l>-^t + 0- (25) 

The vertical velocity Vz is zero in the particular solution (see case B of §2.3); this confirms our previous ad hoc assumption 
(see also Appendix B). As discussed in §2 and it becomes apparent from the equation above, integral-sign warps cannot be 
maintained because of the wrapping action of the disc's rotation. The reason is simple; in a magnetically-driven warp, the line 
of nodes will be described by the locus of those elements for which the vertical magnetic force vanishes, i.e. Bn = 0. Since 
moving with a fluid parcel, Bn is constant in time, the line of nodes will wrap halfway round the galaxy in a time 
7r 

Tnode. = ^^^^^ _ ^^^^^ , (26) 

where Ri and R2 represent the minimum and maximum radii of the warp. Evaluating Eq. (26) for the Milky Way with 

Q{R) = Vo/R, vo — 220 km s~^, we infer Tnodcs ~ 0.32 Gyr, where we adopted Ri ~ 9.5 kpc and R2 ~ 16.5 kpc from Burton 
(1988). Thus, the winding problem in this case is even more severe than in the non-magnetic case for which the wraping time 
of the warp is of the order of ~ 2 Gyr for a quasi-spherical halo (Kahn & Woltjer 1959; Binney 1992). It can be readily shown 
that the line of maximum elevation will also tend to wind up as the line of nodes does, changing the initial assumed warped 
appeareance beyond recognition. The wrapping of the 7i-component of the magnetic field, which is responsible for the winding 
of the warp, is a consequence of a more general result, proved in Appendix B, that no stationary configuration exists for a 
tilted disc in differential rotation under strict frozen-in conditions. 

The increasing tangled magnetic field will act as a pressure in the vertical direction, resulting in an increase in the 
thickness of the disc, and as a tension in the radial direction (e.g., Sanchez-Salcedo & Reyes- Ruiz 2004), which could lead 
to a radial transport of angular momentum by the amplified stress term B^Br (e.g., Kundt 1990). Eventually, because of 
the strong winding of the magnetic lines, magnetic diffusion becomes important (Parker 1973). Hydromagnetic waves in the 
disc will become nonlinear and their energy will serve to heat the disc. This line of arguments strongly suggests that the 
contribution of a torque by the extragalactic magnetic field to a systematic deviation from the mean galactic plane is null at 
the present age of galaxies in the local Universe. 

Not all the galaxies rotate differentially. The observed H I disc of some dwarf galaxies rotate as a solid rigid Q{R) =const 
as occurs in IC 2574 (Blais-Ouellette et al. 2001). For these galaxies, magnetically-generated warps could survive. An obvious 
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result of our analysis is that, even if there was a uniform large-scale intergalactic magnetic field conductive to bend the gas 
discs, a correlation between the orientation of warps as that found in Battaner et al. (1991) is not expected because the line 
of nodes and the lino of maximum height of the galactic warp will rotate with angular velocity fl, different from galaxy to 
galaxy. Therefore, models aimed to study the strength and direction of intergalactic magnetic fields that are based on possible 
correlations between the warps of different galaxies, as those given in e.g. Vallee (1991), are not sustained in any physical 
basis. 

So far, only rotation has been considered. In the next Subsection we will discuss if other effects could alleviate the 
magnetic winding problem. 



3.2 Turbulent diffusion and the a-effect 

Interstellar gas is probably turbulent. On the one hand, turbulent diffusion of magnetic fields can act as a mechanism for 
enhancing transport rates above their molecular values. On the other hand, turbulence together with large-scale flows can 
convert small-scale fields into large-scale fields. The simplest form of the mean field dynamo equation that retains the basic 
physics is 

^ = Vx(t;xB) + Vx {aB - j^tV x B), (27) 

where a w — Z^fi • V ln{pv')F{Bm, f2) and rjt « lv'G{Bm., f2)/3, with I the correlation length of the turbulent motions v' , Bm. 
is the large-scale mean magnetic field, and F and G are certain quenching functions (e.g., Beck et al. 1996, and references 
therein). The a-term represents the effect of cyclonic motions on the mean field. 

Since the disc is warped with vertical displacement h, it is convenient to introduce the variable Z as the vertical distance 
to the midplane of the gas disc, Z = z — h. For simplicity, we assume that the disc is defined by sharp boundaries at Z = ±H/2 
and further introduce the aspect ratio of the disc, A = H/2R. Specifically, we choose H/2 ~ 1 kpc and R = 15 kpc, so that 
A ~ 6.5 X 10"^ 

As the disc is embedded in a constant external field, Be, we expect the same field in the upper and lower surfaces of the 
disc by symmetry: B{Z = H/2) = B{Z = —H/2). Given the even symmetry of the dynamo equations, the components 
and Br must be symmetric with respect to Z: 

B^{Z) = B^{-Z), (28) 
Br{Z) = Br{-Z). (29) 

The ratio between the time-scales of magnetic diffusion across the disc and along the radius is of order of A^. This implies 
that the field distribution across the disc is established relatively rapidly, while the radial distribution evolves on a considerably 
longer time-scale. Consequently, in the absence of the a-efi^ect, B^ obeys a pure diffusion equation dB^/dt ~ r)t^^Bz, where 
the magnetic diffusivity may depend on radius, r]t(R). Diffusion of the plane-parallel ambient magnetic field into the disc 
can counteract the undesirable effect of differential rotation if the appropriate local Reynolds number TZrn = i^H^ /'irjt < 1 
(e.g., Moffatt 1978), implying that the characteristic diffusion time-scale, td{R) = H'^{R)/ir]t{R), at any radius, is smaller 
than local rotation period, i.e. td(R) < l/il{R). In particular, at a radius Ri, td(Ri) ^ 0.1 Gyr. If we turn off the a-effcct, 
turbulent diffusion permits slip of the field in the z direction and tends to make the field more uniform and straight (see 
Fig. 2b). After a few local diffusion time-scales, Bz{R) will take a value equals to the intergalactic value. Hence, the magnetic 
tension responsible for the warp, which is oc dBz/dR, drops monotonically to zero in a few local td- Therefore, turbulent and 
ambipolar diffusion can delay the winding process of magnetically-driven warps but at the price of decreasing the amplitude 
of the warp to a negligible value in a time-scale comparable to Tnodes, i-e. in less than 1 Gyr; the inclusion of any magnetic 
diffusion aggravates the problem of the magnetic support of galactic warps. 

To show that the a-term is not suitable to change the dramatic damping effect of turbulent diffusion, consider the mean 
vertical value of Bz in the disc 

Fz{R,4>) = T7 / BzdZ. (30) 

^ J-H/2 

From the ^-component of Equation (27) 



dt H J_H \R dR ' R 

4 f4 + c.B,dZ-4^l aBRdZ + ,^V^F ' 



H \R dR _H HRdS _H " " ' H dZ 



(31) 
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Because and are even in Z (see Eqs. 28-29) and a is odd, the integrals involving the product a and or the product 
a and Br vanish and, hence, 



dt H dZ 



(32) 



Equation (32) indicates that the diffusion of the component Bz across the cap surfaces of the disc aX Z = ±H/2 follows 
unimpeded even in the presence of the a-effect. 

In conclusion, in order to avoid a fast reduction of the amplitude of the warp to unobservable values, the magnetic 
diffusivity should be significantly quenched relative to the canonical turbulent diffusion value, rjt. 



3.3 The magnetic field outside the disc 

Regular magnetic fields have been detected not only in the discs of spiral galaxies but also far above them, in the surrounding 
gaseous medium. The reason is that both the magnetoionic disc, which is thicker than the H i disc, and the gaseous halo are 
magnetically active. Because we were only interested in the magnetic term responsible for the warp formation and since warps 
are observed in the very cold Hi disc, the details of the evolution of the complete magnetic field outside the disc are not 
important. A thorough analysis of the magnetic field structure is out of the scope of this paper. For completeness, however, 
we will comment on the magnetic field outside the disc. 

Initially the magnetic field was supposed to be uniform and ccjual to Be, the intergalactic field. The field trapped by the 
galaxy will undergo the action of chaotic motions and differential rotation. All we can say, at this point, is that the outer 
boundary condition should be that the field Be at infinity is undisturbed, i.e. B — > Be as r — > oo. 

In order to study the field structure in the magnetoionic disc, we need to implement boundary coiiditioiis at the surface 
of the disc in Eq. (27). The teneous and partially ionized gas in the corona suggests force- free conditions for the medium 
outside the disc. Observations suggest that this is not the case. Coronal gas is in chaotic motion, and acquires a nonzero mean 
helicity due to galactic rotation (e.g., Sokoloff & Shukurov 1990). Since the dynamo can be active not only within the disc but 
also in the gaseous halo, one should resolve, in principle, the disc and halo fields in a self-consistent calculation. However, the 
magnetic diffusivity in the halo is ~ 50 times that in the disc, which indicates that it is meaningful to consider the magnetic 
evolution in the galactic disc by neglecting the dynamo in the halo and adopting vacuum boundary conditions (e.g., Poezd 
et al. 1993; Bardou et al. 2001). Numerical simulations by Moss & Brandenburg (1992) confirm that the vacuum boundary 
conditions can be applicable in this case. In this context, a vacuum is understood as a medium in which electromagnetic 
induction effects are significantly weaker than the magnetic diffusion. 

If the disc and the corona are rotating as a rigid body, the evolution of the regular field in the disc would be quite simple; 
q; ~ and diffusion effects lead to the external field Be to penetrate into the disc; the ultimate steady-state is a uniform field 
and equal to Be- Differential rotation makes the situation more complex. The response and development of the field lines 
due to rotation is a classical problem in MHD. The full time-dependent problem of a rotating cylinder embedded in a perfect 
insulator weis solved by Parker (1966). The interesting reader is referred to this paper for a visualization of the lines of force 
and its different behaviour with the Reynolds number (see also Moffatt 1978, Section 3.8). The evolution of the field lines in 
the standard cooling flow scenario for the formation of disc galaxies can be found in Rognvaldsson (1999). 



3.4 Other caveats 

In this discussion we have assumed that the center of mass of the galaxy is at rest in the initially plane-parallel intergalactic 
magnetic field, as a prerequisite for galaxies to be integral-sign warped in the magnetic scenario. For galaxies moving through 
the intergalactic medium, shear terms will distort the field lines in the outer parts of the disc and, therefore, the field lines 
cannot be considered plane-parallel any longer, but a magnetic barrier of increased field strength would exert a pressure on the 
disc (e.g. Konz ct al. 2002), taking the classical form of a rim. In the next Section we will discuss whether the hydrodynamical 
force due to the accretion of intergalactic gas, which is possibly more important than the magnetic forces in moving galaxies, 
can contribute to the formation of warps. 



3.5 Comparison with previous work 

Battaner et al. (1990) proposed the galactic warps may have a magnetic origin since there exists a magnetic field distribution 
liable to produce a m = 1 torque in a slab of slightly ionized gas. Interestingly, the magnetic configuration presents lines that 
open to the space and connects at infinity with a uniform magnetic field, which was identified as the ubiquitous intergalactic 
magnetic field. A detailed calculation of the magnetic field structure required to reproduce the amplitude and shape of the 
warp in NGC5055 was given in Battaner & Jimenez- Vicente (1998). Their analysis was carried out by considering the Lorentz 
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force in the vertical equation of motion and the condition V • i? = 0, although ignoring the induction equation for the field, 
or more precisely, the field was assumed to be in a steady-state configuration, i.e. dB/dt = 0. Under these assumptions, the 
predictions of this model arc: (1) the postulated intergalactic field is of order of a few fiG, (2) a spatial correlation between 
the direction of the warps in neighbouring galaxies and (3) a shift between the stellar and gaseous warps. 

However, it is clear that the magnetic field in the disc will be advected with the fiow. We have studied this process by 
including the electromagnetic induction effects in the postulated field. The assumptions of Battaner ot al. (1990) and Battaner 
& Jimenez- Vicente (1998) are only valid for slightly turbulent systems with a rigid rotation, such as the magnetic configuration 
is not destroyed by differential rotation or turbulent diffusion. Even for galaxies rotating as a solid body, however, the rotation 
of the galactic discs excludes any alignment between the orientation of warps of nearby galaxies by a large-scale intergalactic 
field. For discs in differential rotation, we show that as the magnetic field evolves due to the stretching term, warps become 
twisted and destroyed in less than two or three rotation periods. The increasing tangled magnetic field will act as a pressure, 
forming a unwarped thick H I disc at later times. The thickness of the H I disc will depend on the turbulent diffusion time-scale. 
If it is shorter than the time-scale of winding, the warp damps its amplitude before the disc gets thicker. In our analysis we 
have adopted a prescribed large-scale velocity field ignoring the magnetic back reaction of the field onto the gas dynamics, 
which could load to a radial inflow of gas. 

If galaxies travel with a peculiar velocity through a medium permeated by an ordered magnetic field of ~ 1 jiG, the 
magnetic barrier of increased field strength would contribute to the ram pressure in forming a U-warp. Therefore, we suggest 
that features in the morphology of the outer H I disc can shed light on the strength of the regular intergalactic magnetic field. 



4 GENERATION OF WARPS BY AN INTERGALACTIC ACCRETION FLOW 
4.1 Revisiting the accretion model: motivations and caveats 

It is well-known that the dynamical pressure acting on a galaxy travelling through the intergalactic medium (IGM) could 
bend the outer parts of spiral galax;ies, especially in cluster galaxies for which the interstellar gas can be even stripped from 
the galaxy (Gunn & Gott 1972). Several galaxies in clusters are reported in the literature whose distorted and asymmetric 
extraplanar H I gas is claimed to be suggestive of ram pressure stripping by the intraclustcr gas, rather than accretion. The 
assumptions behind the steady-state accretion model of small vertical perturbations, being the disc a perfect absorber of the 
incident intergalactic gas, do not apply for cluster galaxies moving supersonically in a teneous hot medium. The aim of this 
Section is to test the accretion model as an explanation of the origin of warps in relatively isolated galaxies or galaxies in 
groups. Therefore we restrict our analysis to these systems since the approximations and results derived in this Section cannot 
be extrapolated to study the extraplanar gas components observed in cluster galaxies. 

It is illustrative to estimate the effect of the ram pressure of a wind that collides directly onto the disc. In order for ram 
pressure to bend the gas disc a height h, with h <^ R, above the mean plane, the following order-of-magnitude condition has 
to be met: 

Poavt, > SgZ/rf/l, (33) 

(e.g., Gunn & Gott 1972; Marcolini et al. 2003) where poo represents the mean gas density of the IGM, Voc the velocity of the 

galaxy through the IGM, Eg the disc's gas surface density, 7? the galactocentric radius, i^t the disc vertical frequency and Viot 
the circular speed. For an exponential disc, Vd{R) — ^{R)^/^ + 22{Rd/R)''^ (Binney 1992). For the following reference values: 
Vrot ^ 200 km s~^, Voo = 200 km s^^ (see, e.g.. Bureau & Carignan 2002), i? = 15 kpc and Eg = 1 M© pc~^, we require 
poo > 10~* cm~^ to have a displacement oih ~ 500 pc. The present-day mean density of the IGM beyond galaxy haloes in the 
Local Group and in small groups of galaxies is thought to be < 10~^ cm~^ (e.g., Murali 2000; Bureau & Carignan 2002; Grebel 
et al. 2003). With this intergalactic volume density only /i < 50 pc could be expected for the above reference values, which 
seems to be too low to explain warps in large spiral galaxies. As mentioned at the beginning of this Section, since the values 
of Poo and Voo are much larger in clusters of galaxies, we cannot exclude the possibility that extraplanar H i and warps seen 
in some galaxies showing signatures of ram pressure stripping are induced by the interaction with the intracluster medium. 
Good candidates are the galaxies NGC 4548 in the Virgo cluster (VoUmer et al. 1999) and HoII in the M81 group (Bureau & 
Carignan 2002; but see also Bureau et al. 2003). To our knowledge, numerical experiments of ram pressure stripping do not 
show any evidence of the formation of integral-sign warps (e.g., Quilis, Moore & Bower 2000; Marcolini et al. 2003; Roediger 
& Hensler 2004). Other possibilities for the explanation of distorted velocity fields in cluster galaxies have been suggested 
(e.g., Vollmer et al. 1999). 

In the derivation of Equation (33) it was implicitly assumed that the flow behaves as a continuous gas medium. Inter- 
estingly, LC show that if the flow is a gas made of discrete coUisionless undisruptive clouds, the gravitational focusing by 
the galactic potential would enhance the density of the flow by an important factor when it impacts the disc. In fact, they 
obtained that the northern warp of the Milky Way {h ^ 2 kpc) can be explained by the infall of matter through a subsonic 
wind of 100 km s~^ with a density of ~ 6 x 10~® cm~^. 
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Figure 3. Graphical representation of a vertical cross section at 2: = showing the warped disc. Panel (a): Precession of the galactic disc 
with frequency lu around v^c as a single unit (model LC). If the disc is embedded in the flattened potential of a massive dark halo, which 
we take to be aligned with the coordinate system (panel b), the direction of the incident flow is not any longer the axis of symmetry of 
the system. Both the dark halo and the central disc are taken as frozen and unresponsive bodies because are so massive that they are 
very reluctant to prcccss. Contrary to panel (a), the outer rings cannot achieve large vertical heights with respect to the equatorial plane 
of the dark halo (z = 0) because the dark halo and the inner disc try to twist the disc back into alignment with the coordinate system. 
The steady-state configuration corresponds to a; = (see text). 



For a flow approximately parallel to the rotation axis of the disc (face-on wind), the shape of the outer gas disc should 
resemble an a^xisymmetric U-rim instead of an integral-sign warp. If the flow is oblique, the impulse transfer has an azimuthal 
dependence and, hence, the vertical distortion h is a combination of the m = mode, giving a U-shape profile, and the m = 1 
mode (S-shaped proflle). The success of the ring-tilted decomposition proposed by Rogstad ot al. (1974) in reproducing the 
observed morphology and kinematics of real warps suggests that the m = 1 component is indeed dominant with respect to the 
m = one. Hunter & Toomre (1969) noticed that in the wind model by Kahn & Woltjer (1959), warps should be asymmetric, 
with the maximum gas displacements on the two sides of the Galaxy in the ratio 3 : 2. The similarities between the wind 
model and the accretion model suggest that the warps in the accretion model should be also cisymmetric. Unexpectedly, LC 
found that the S-amplitude (m = 1 mode) is much greater than the U-amplitude (m = mode) except when the flow vector 
is within 5° of the spin axis of the galaxy. They claim that if all galaxies had Milky Way parameters around 1 in 1000 should 
have a U-shaped warp (Beckman et al. 2003). This result is very counterintuitive indeed. 

In the stationary models of LC^, the galactic disc as a whole, and not only the warp, precesses as a rigid unit with a 
no-null angular velocity cj around Woo, the axis of symmetry of the system (sec §2.1 and Fig. 3a). As discussed in §2.1, since 
the disc is rather massive and rotates fast, u) must be small but different from zero. Assuming that the dark halo is perfectly 
spherical, the precessing frequency about Voo, is given by = Tacc/C-'^disc cos^ ^0), where Tacc is the external torque acting 
on the disc caused by the accretion flow of intergalactic gas and cos So = sin(uoo, idisc). After estimating Tacc following LC, 
one gets that the precession rate around Voo for an angle of 45° is of ~ 0.18 km s~^ kpc~^ (see also LC). The dynamical 
correction for this effect is apt to be ignored and hence we can treat the inner disc as a fixed external potential. Under that 
approximation and considering only the evolution of small vertical displacements from the midplane z = defined by the 
fixed central disc, the torque applied by the fiow will be constant in direction and magnitude. This situation corresponds to 
case A of §2.3 and, therefore, a; = is expected in the stationary configuration. Similarly, if we include a massive dark halo 
flattened along the rotation axis and treat it as a time-independent background potential, the outer rings of the disc cannot 
present a steady-state precession around Voa, as occurs in a spherical halo (see the disc at t = n/uj in Fig. 3a), because the 
halo twists the disc back forcing it to stay close to the halo midplane z = (see Fig. 3b). Once the halo and central disc 
arc treated as time-independent external potentials, we will confirm in §4.3 our expectation that a; = in the steady-state 
configuration. However, this result has a simple explanation, as follows. 

The change of the strength of the angular momentum of each ring, £ring, is 

— 2Tring ' iring, (34) 

where Tring is the total (gravitational plus nongravitational) torque. In a stationary configuration dL^^^^/dt = and hence 
either Tring = (and hence the rings do not precess) or the torque is always perpendicular to iring- As already discussed, for 
the time-scales of interest, the torque applied by the flow can be considered constant in direction and magnitude, whereas the 



^ We must distinguish here between the halo of coronal gas and the halo made of coUisionless dark matter particles. In the accretion 
model, the coronal gas is ignored whereas the dark halo may be present. The interaction of the accretion flow with the dark halo is only 
through the gravitational force and, hence, the flow can rea<;h the disc. 
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gravitational torques are always perpendicular to iring- Therefore the condition Xring always perpendicular to iring cannot 
be met and, in fact, the only stationary configuration corresponds to the particular solution of case A in §2. In this particular 
configuration there is a balance in each ring between the non-gravitational torque produced by the accreting fiow and the 
external plus internal gravitational torques (see Eq. (40) and Appendix A) . 

The halo not only alters the way in which the rings precess. It has been recognized for a long time that the flattening of 
the dark halo is a crucial matter in order to understand the evolution and shape of vertical perturbations in discs, because 
the amplitude of the warp depends on the contribution of the dark halo to the vertical restoring force, keeping the rings more 
tightly bound in a flattened halo. 

In the next Sections we derive estimates of the amplitude, warp asymmetry and evolution of the tilted disc using a 
formalism that allow us to treat all the m-modes in a self-consistent way, taking into account the halo flattening, which 
is essential in modelling the vertical distribution. Our procedure is completely different to the approach followed by LC, 
motivated to take a different look at the problem. The mode decomposition will help to better understand the underlying 
physics of the situation. For instance, our calculations can tell us when the approximation of solid rings is correct. In the 
following Subsections we will discuss strictly the predictions of our calculations to test the viability of the accretion scenario. 

4.2 Assumptions 

The calculations start from the following assumptions: 

(1) The accretion of an intergalactic flow directly onto the disc is the main physical agent responsible for the formation 
of warps. 

(2) The intergalactic flow is composed of purely baryonic matter and contains a negligible amount of collisionless dark 
matter. 

(3) The flow is made of discrete collisionless clouds with a null velocity dispersion at inflnity. 

(4) The fiow is plane-parallel at infinity, with the direction and velocity constant in time. 

(5) The clouds move in hyperbolic orbits, unaffected by the drag with the hot coronal gas halo, in a galactic potential 
described by a point-mass. For the Galaxy the adopted mass is 2 x 10^^ M© , which is approximately the dynamical mass of 
the Galaxy within a sphere of 16 kpc in radius. All the mass that intersects the disc is perfectly accreted. 

(6) The inner disc and halo are frozen or non-responsive, i.e. they are modeled as time-independent external potentials; 
we do not include the response of the inner disc and halo to the tilt of the disc. The central disc and the halo are aligned. 

(7) The amplitude of the warp is so small that the force caused by the intergalactic flow does not depend on it. This 
assumption together with assumption (6) imply /acc(-R, ^o), where 9o is the angle of incidence of the flow. Consequently, 
given ^0, the collisional torque exerted by the flow is constant in magnitude and direction. 

(8) The gravitational torque exerted by the flow into the disc is neglected. We consider only the torque produced by the 
collision of the clouds with the disc. 

(9) The stellar and gas discs are gravitationally tight enough that the differences between stellar and gas warps are small. 
All these assumptions wore also adopted by LC except assumptions (6) and (7). We have estimated that the neglect of 

the precession of the entire disc to be a reasonable approximation (see §§2.1 and 4.1). Assumption (2) is somehow related 
with assumption (8); gravitational torques by the accretion flow are expected to be small compared to the collisional torque 
for the parameters under consideration. Concerning assumption (3), a null velocity dispersion of the clouds is the simplest 
assumption. From basic grounds, it can be seen that the relative amplitude of the m = 1 mode respect to the m = mode 
depends on the velocity dispersion of the clouds. Regarding assumption (5), it is likely that clouds disrupt before reaching the 
inner disc. For these clouds the hydrodynamical drag with the gaseous corona cannot be neglected (e.g., Benjamin & Danly 
1997). However, it is equally possible that certain dense clouds can achieve the far-outer parts of the disc before disruption, 
for a plausible range of parameters of the clouds. 

4.3 The accretion model: Equation of motion and the mass model 

The vertical dynamics of a parcel of material of the disc with vortical displacement h[R,(j),t) and surface density S(i?), 
embedded in a dark halo with gravitational potential $h, is given by the equation of motion 



where Q,{R) is the circular frequency at radius Ji, = d'^^-h/dz^\z=o is the vertical frequency of the dark halo, while 
/self (i?, 0, t) and /acc(-R, 0) are the vertical forces per unit of mass on the parcel of gas due to the self-gravity of the distorted 
disc (/self) and due to the intergalactic accretion flow (/acc). 

The dark halo is assumed to be frozen, i.e. it behaves as a flxed potential, and flattened along the symmetry axis of the 
disc with axis ratio q, typically q w 0.75, with the classical pseudo-isothermal density proflle: 




(35) 
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Table 1. Parameters for six models (columns 2-3), together with the amplitude ratio between the two sides of the warp at R = 22rc 
(column 4) and the shape of the warp (column 5). The two latters were derived by solving Eqs. (41) and (42). S denotes integral-sign 
warps and L stands for one-sided warps (see Fig. 4). Type I refers to warps with increasing amplitudes towards the edge of the disc. 
Asterisk indicates that a velocity of Voo = 300 km instead of 100 km was used. 



Model 


? 


eo(deg) 


7l 


Shape 








at R = 22rc 




I 


0.75 


45 


1.54 


S, type I 


II 


0.9 


45 


1.39 


S, type I 


III 


0.75 


30 


1.37 


S, type I 


IV 


0.75 


60 


1.97 


S, type I 


V 


0.75 


75 


5.43 


L, type I 


VI* 


0.75 


45 


2.02 


S, type I 



Ph(i?,a) = <^ ^ ^, . (36) 

I otherwise, 

with = B? + /q'^ , the core radius of the dark halo, rt = SOtc and Too = 65'rc, those last values were selected just to 
facilitate comparison with previous models, as in Binney et al. (1998), but they do no play any noticeable relevance in the 
accretion model provided they are larger than the maximum radius of the warp R2 (see assumption (5)). 

The warps are observed in both the stellar disc and the gas layer and, where they can be compared each other, the 
warped planes are almost the same (Casertano et al. 1987; Cox et al. 1996). Hence, within approximation (9), the surface 
density of the disc should include both stars plus gas. We adopted a truncated exponential profile for the surface density with 
scalelength R^, i.e. E(i?) = So cxp(— i?/i?d), whenever R < Rt, with Rt the truncation radius. 

The vertical force /self weis given in Eq. (5). Following Binney et al. (1998), we have adopted a softening radius zo = O.IR^. 

The physical basis for the warp formation in the model of intergalactic matter accretion relies on the fact that /acc may 
exert a net torque in the disc because the flow is redirected and gravitationally focused by the potential of the galaxy. With 
no loss of generality, we will assume that the accretion flow comes from below the equatorial plane, i.e. Vz > 0, with a velocity 
at infinite distance of Voo = {vx,0, Vz) in the galaxy rest frame. We define the angle ^0 between the direction of the incoming 
gas flow at infinity and the reference plane, i.e. sin6'o ~ Vz/vao- A flow parallel to the spin axis of the galaxy has ^0 = 7r/2. 
/acc was derived by LC in the approximation (5). The formula for /acc is given by 

/acc(i?,0)=^°°''' 



S(i?) 



where b{R,(j)) is the impact parameter and (pbi4>) is the polar angle in the plane perpedicular to Vao- The expressions for 
b{R,<f)) and 4>b were given in LC [their Eqs (38) and (34), respectively]. 

The imposed force /acc can be expressed as a sum of azimuthal Fourier terms of different amplitudes, 

00 

/acc(Ji,0)= ^/r.(i?) cos m(j). (38) 

The steady response of the disc will be also a sum of Fourier modes with amplitudes ho{R), hi{R), h2{R) and so on. The rings 
can be treated as rigid provided that the amplitude of the displacements for components m = and m = 1 exceed those for all 
other wavenumbers. We have checked that this is the case except for small values of 60. In fact, for 60 < 15°, the /i-amplitudes 
of m = 2 mode is comparable to that for the m = component. Hence, for 9o > 15° , we take /acc {R, 4>) — fo (R) + fi (R) cos 4> 
and omit the analyses for other m-components, those analyses being analogous. 
The steady response of the disc has the form 

h{R,ct>,t) = ho{R) + hi{R)cos{u;t-ct> + ct>o{R)), (39) 

where u is the pattern speed about the 2;-axis, ho{R) and hi{R) arc the amplitudes for the m = and m = 1 modes, 
respectively, and (f>i,){R) a phase function which describes the twist of the line of maximum heights of the pure m = 1 mode. 
The observed integral-sign warps, or S-shape warps, correspond to the second term (m — 1), whereas the first m = mode 
is responsible for the generation of the cup-shaped symmetry between the two sides of the warp. 
Substituting the above expressions for h, /self and /acc into equation (35) yields: 



MR) = 0, oj = 0, 



(40) 
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fo{R) = vlhoiR) + Gho{R) / T,{R')H{R,R')R'dR' 




-G i Y:{R')H{R,R')ho{R')R'dR', 



(41) 



fRt 

fi{R) = {i^l - Q')hi{R) + Ghi{R) / T.{R')H{R,R')R'dR' 




-G / J:{R')J{R,R')hi{R')R'dR' , 



(42) 



with H{R,R') and J{R,R') as defined in Eqs. (15) and (16), respectively. 

As wc anticipated in §4.1, the zero pattern speed a; = 0, which means that the rings does not precess around the z-axis 
in the steady forced evolution, is a consequence of our assumptions (6) and (7). Under these approximations the form of the 
imposed time-independent forcing is /acc(J?, </>) = fo{R) + /i(J?)cos0, whose m = 1 mode corresponds to case A of §2 (see 
also Appendix A). For the same reason, the twist function vanishes, </>() — 0, implying that the line of maximum heights is 
straight; the direction of the torque produced by the accreting flow provides a direction along which the warp's line of node 
may naturally align. Unfortunately, the predicted value w = is not in accord with the pattern speed inferred from the stellar 
velocity distribution in the warped Galactic disc of a; < —20 km s^^kpc^'^ (Smart ct al. 1998). 

The integral equations (41) and (42) can be solved by treating the continuous disc as a collection of massive rings and 
using standard methods. We will use units such as G = rc = 1 and the halo mass is unity. For the disc mass we choose 0.1. 
Note that our mass model is identical to that adopted by Binney ct al. (1998). The disc is made up by 60 rings uniformly 
spaced in radius between R = 9rc and R = Rt = 22rc. The results for four representative models with different angles of the 
wind and flattenings of the dark halo, as described in Table 1, will be explored and used as a test for the accretion model in 
the next Subsection. 



4.4 Amplitude and shape of the warps: Results, interpretation and discussion 

It is convenient to start our discussion by considering the baryon density and flow velocity required to explain the observed 
amplitudes of galactic warps. Later on, we will consider the observable asymmetries thereby detectable in H I velocity contour 
maps and in edge-on galaxies, as a result of the coexistence of the modes m = and m = 1. A comparison of predictions with 
observations will be used to validate or rule out the accretion model. 



4-4-^ Amplitude 

Figure 4a displays the height of the disc respect to the symmetry plane of the halo versus the radius of the disc along a 

cross-section cut through the warp for some models specified in Table 1. In all the models, we have used the same parameters 
for the flow than in LC, Woo = 100 km s^^ and poo = 6 x 10~^ cm~^. This density is likely a generous assumption (see 
§4.1). For a galaxy like the Milky Way, Rd = 4.3rc = 3 kpc. It turns out that the observed Galactic warp amplitude can be 
reproduced with the adopted parameters. We confirm the results of LC derived in a very different framework. 

For cluster galaxies, both poo and Va^ can be significantly larger than the values quoted. Nevertheless the predictions of 
the accretion scenario are not longer valid for these galaxies because the assumptions behind the accretion model are not a fair 
representation of the coUisional interaction between the intraclustcr medium and the interstellar gas (e.g., Sanchez-Salcedo 
2004). Numerical calculations of the ram pressure stripping process are required to explain the distorted velocity field and 
extraplanax (usually resembling a U-rim) distribution of the gas observed in these galaxies (e.g., Vollmer et al. 2000). Since 
the problem of the persistence of warps is circunscribed to relatively isolated galaxies, we will omit any discussion regarding 
the origin of warped galaxies in clusters. 

Fixed Poo and Voo at the preferred values, a scaling relation between the size of the galaxy, Ri, defined as the last Hi 
measured point and warp's amplitude could be expected. The accretion force given in Eq. (37) depends on the parameter 
GM/Rv^. To a good approximation and for the range of interest, it can be fitted by /acc oc {GM / Rvi^Y "^ ~ (vrot/voo)^'^, 
where we have used that M is the dynamical mass within Ri, as adopted in previous Sections. The displacement h at the 
end of the warp hi scales according to /acc ~ hiv'^^tR^^ (see Eq. 33). Therefore hi ~ v'^^^Rf and we would expect a positive 
correlation between the amplitude of the warp and two parameters; galaxy radius and galaxy mass. In contrast, Castro- 
Rodriguez et al. (2002) found a slight anticorrelation between the observed amplitude and these parameters. In fact, the 
amplitude of warps in dwarf galaxies cannot be accounted by the accretion model itself. 
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Figure 4. Shape of the warp along a cut through the (straight) line of maximum heights for the different models of Table 1 (panel a) 
The radial dependence of 7, which is a measure of the warp amplitude asymmetry is shown in panel b. For making estimates in the 
context of the Galaxy, rc — 0.7 kpc. 

4-4-2 Shape and morphological asymmetries 

As it becomes apparent in panel (a) of Fig. 4, the accretion model predicts that all the S-warps must be of type I in the 
terminology of Sparke & Casertano (1988), since the inclination of the outer warp increases from the center to the edge of 
the disc. This result is mainly a consequence of the dependence of the force term with the surface density /acc oc (see 
Eq. (37)). If the flow is smooth and accretion is the only cause of warps, it is very difficult to explain why it is very common 
to find that one side turns back up toward the inner-disc plane after reaching its maximum vertical height (Burton 1988; 
Garcia- Ruiz et al. 2002), implying that the amplitude of the m = component relative to that of the m = 1 mode increases 
with galactocentric radius. As an example, the curve of the asymmetric warp of the Milky Way is shown in Fig. 5 from Burton 
(1988) with Rq = 8 kpc. For our Galaxy, the ratio between the amplitude of the mode m = over the amplitude of mode 
m = 1 increases from ~ 0.1 at 14 kpc to 0.7 at 18 kpc. This radial behaviour is unlikely in the accretion model as the m = 1 
response becomes dominant at the outer parts due to the near equality of the orbital and the vertical oscillation periods at 
successively larger radii. As it is illustrated in Fig. 4a, the warps never curve back to the disc-plane. 

In the accretion model warps should be asymmetric (see panel (a) of Fig. 4). Given the mass model of a certain galaxy, 
the asymmetry of the warp depends on 9o, Vao, q, and the velocity dispersion of the clouds, which was assumed to be zero. 
Smaller the angle 9o is, more antisymmetric (i.e., they are closer to a integral-sign shape) the warp becomes, but even for 
9o = 30°, the amplitudes between the two sides in the warp are significantly different (see Table 1). It is interesting to note 
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Figure 5. Shape of the warp of the Milky Way taJjen from Burton (1988) with Rq = 8 kpc. 

that for the model with = 45°, the ratio between the amplitudes fi/fo ~ 1, fairly independent of R in the range of interest, 
but the disc is more responsive to the m = 1 mode, resulting in hi /ho « 5 at 22rc. The reason is that the m = amplitude 
ho oc u^^ (in the notation of §2), whereas hi oc {u^ — In panel (b) of Fig. 4 we have plotted the radial dependence of 

the value 7(-R). defined as the ratio between the amplitudes at the same 7?, as measured from the central galaxy plane to the 
peak of Hi column density, on both sides of the warp with the convention that 7(i?) > 1. For model IV, 7 — > oo at a certain 
radius near R ~ 16rc. This divergence is because at that radius the warp crosses the equatorial plane of the inner disc at 
one side. For angles 9o > 75° warps bent with a L-shape (see model V in Table 1) or even a bowl-shape. Hence, we expect a 
percentage of occurrence for these types of warps of ~ 3.5%. This fraction may be larger if the clouds have no-null velocity 
dispersion at infinity. In fact, if the flow velocity has its origin in the motion of the galaxies, one expects similar peculiar 
velocities for the clouds. In that case, i.e. clouds having an intrinsic velocity comparable to Voo, the percentage of bowl-shaped 
warps increases by a factor ~ 2. 

Observations show that the integral-sign warps are not perfectly symmetric (Sancisi 1976; Burton 1988; Bosma 1991; 
Garcia-Ruiz et al. 2002). The natural question that arises is whether the accretion model can explain the degree of asymmetry 
of warps of type I. To check this possibility, we have selected from the literature a sample of 16 galaxies with type I warps 
and estimated ji, the value of 7 at the end of the warp. We have not included those galaxies that are strongly perturbed by 
the interaction with nearby neighbours or with very small warps. The selected galaxies are: M33 (Rcakes & Newton 1978; 
Corbelh & Schneider 1997), NGC 4244 (Sancisi 1976), NGC 300 (Rogstad et al. 1979), NGC 2841 (Begeman 1987), NGC 4013 
(Bottema 1996), UGC 7170 (Cox et al. 1996), ESQ 123-G23 (Gentile et al. 2002), NGC 5055 (Battaglia et al. 2005), U6964, 
U7774 and U8711 from Garci'a-Ruiz et al. (2002) and NGC 2685, NGC 3718, NGC 5204, UGC 3580 and NGC 2541 (Jozsa et 
al. 2004 and private communication). All the galaxies of this sample present integral-sign shaped warps with 7; < 2.15, and 
most of them are so symmetric that no difference between the amplitude of both sides has been reported. In order to increase 
statistical significance, we have considered the distribution of the 7;-values for the sample of 137 warped galaxies observed in 
the optical band (Castro- Rodriguez et al. 2002). We find that only 36% of these galaxies present optical warps with 7; > 1.25. 
Note that since stars extend less than the gas layer, the stellar warp is expected to be even more asymmetrical (see Fig. 4b) . 
The existence of galaxies with small values of 7; is difficult to explain if accretion is the main and only responsible for warps. 

From studies of the asymmetric shape of the warps, we conclude that accretion flows may play a role in warping the 
discs, but they cannot be the only cause of warps. 

4.5 The vertical distribution in H I volume density 

Sanchez-Salcedo (2004) pointed out that if the force responsible for the warp has a hydrodynamical origin, rather than 

gravitational as in the cosmic infall scenario, the vertically shifted disc must be supported by a pressure gradient. The reason 
is that one can impose the condition of hydrostatic equilibrium as the sound crossing time over the height of the disc at 
the sound speed of 10 km s~^, is much shorter than the period for the gas to cross two successive crests of the warp. 
Consequently, the volume density of gas in a cut along the vertical axis z, at a given (JZ, 0), must present a very asymmetric 
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vertical distribution, but this seems to be in conflict with the observed fairly symmetric distribution of the Galactic warp by 
Diplas & Savage (1991). In addition, a high rate of energy deposition in the cold disc by the ram pressure might produce 
other observable consequences, such as the ionization of the Hi layer (Kenney & Koopmann 1999). 



4.6 The thickness of the outer H I disc 

The scaleheight of the outer disc may provide additional constraints on the strength of the ram pressure acting on discs. 
In addition to the asymmetric distribution in a cut along the vertical z-axis discussed in the previous Section, the external 
dynamical pressure caused by falling clouds, will produce a vertical compression of the Hi disc. This compression may have 
important observational consequences for the outer thickness of the H i disc, such as a reduction and a large-scale asymmetry 
in the Hi scaleheight. The asymmetry is expected as the result of the azimuthal dependence of this external pressure, which 
ultimately is responsible for the S-warp. The recent work by Narayan et al. (2005) implies there is little room for asymmetric 
infall of material at the outer warp. In fact, we suggest that studies of the Hi scaleheight can place strong constraints on the 
accretion rate of intergalactic gas. To illustrate this point, we are going to estimate the thickness of the far-outer regions of 
the Galactic disc for the preferred parameters of the accretion model: 6o = 45°, Poo = 6 x 10~® cm~^ and Voo = 100 km . 
When comparing with observations we will adopt a value Rq = 8 kpc. 

The external pressure due to accretion can be written as a(i?, 0)pooV^, where the dimensionless function a{R,(p) can 
be derived from Eq. (37). At a galactocentric distance of 22.5 kpc a « 40 at = 0, and declines to a value ~ (40/6) at the 
opposite side 4> = n. An upper limit on the thickness of the disc, as measured as the height containing 63% of the Hi surface 
density, based on the equation from equilibrium is: 



L < min 



2c: 



■g ^a'-g 



_7rG(E5 + (c9/c*)E*) ' apsofic. J ' ^ "^^ 

where Eg and E« are the surface density of gas and stars, respectively, Cg is the effective gaseous vertical velocity dispersion, 
including cosmic rays and magnetic fields, and c* the stellar velocity dispersion. The first expression inside the parenthesis 
corresponds to the thickness of an isothermal self-gravitating sheet in the lack of any ambient pressure (e.g., Elmegroon 2002), 
while the second expression gives the characteristic thickness of an isothermal gas layer that is pushed in a cylinder from one 
side by a piston moving with constant acceleration, neglecting self-gravity. The action of the piston mimics the pressure of 
a rain of clouds falling on one cap of the disc. Since both self-gravity and dynamical pressure occur simultaneously in the 
accretion model, Eq. (43) necessarily provides a boundary for L. To evaluate L, we need the Hi velocity dispersion, which is 
observed to be almost constant with radius and is about 9 ± 1 km s~^ in the inner Galaxy (Spitzer 1978; Malhotra 1995). 
Beyond the solar circle, however, the dispersion is not yet measured. Studies of external galaxies show that the observed 
velocity dispersion has a narrow range (e.g., Lewis 1984; Shostak & van der Kruit 1984; Dickey et al. 1990). Aiming at 
conservative estiamtes, we will take a dispersion of 10 km s~^. If Cg ~ 7 km s~^ without contribution from magnetic effects 
(as measured directly by the H l velocity dispersion in some galaxies) , the scalcheights would bo about a factor 2 smaller. For 
Eg = 1 Mq pc"^ and E* « at = 22.5 kpc (Wouterloot et al. 1990) and Cg = 10 km s"\ we get L < 180 pc at = and 
L < 1100 pc at = TT, and half of that for Cg = 7 km s~^. This variation of at least a factor of ~ 6 in azimuth, is not supported 
by observations that show very little net asymmetry in the derived thickness of the Hi layer (Henderson et al. 1982; Merrifield 
1992). Moreover, a thickness of 90-180 pc at ~ 22.5 kpc is in serious conflict with the observed measurements of 2200 pc in 
the Galaxy (Wouterloot et al. 1990) and of ~ 1 kpc in external galaxies (e.g., Sancisi & Allen 1979; Brinks & Burton 1984; 
Oiling 1996). We find that it is not possible to bend the outer disc a kpc- height without making the scaleheight too small. By 
imposing that the outer disc is not appreciably compressed by the cloud rain, the height h derived in §4.4 should be reduced 
by a factor ~ 10. 



5 CONCLUSIONS 

The failed of the theory of bending modes of self-gravitating discs in explaining the persistence of galactic warps (Nelson 

6 Tremaine 1995; Binney et al. 1998) has led to invoke external torques caused by the interaction of the galaxy discs with 
their environment. The intergalactic magnetic field and the accretion of a teneous intergalactic gas onto the disc have been 
suggested as possible causes for galaxies to be warped. The fact that galaxies in poor environments axe even more frequently 
warped than galaxies in dense enviromcnts (Garcia-Ruiz et al. 2002) and other empirical correlations (Battaner et al. 1991; 
Castro-Rodn'guez et al. 2002; Sanchez-Saavedra et al. 2003) make the above formation mechanisms very attractive and it 
seems appropriate to seek ways in which we can test their viability. 

We have shown that magnetically induced warps suffer from an acute winding problem as differential rotation is unimpeded 
under frozen field conditions and hence warps cannot persist. Curiously, the time-scale for winding of the warp is shorter than 
in the non- magnetic case studied by Kahn & Woltjer (1959). It remains unclear why the different warps of neighbouring 
galaxies should be aligned if the magnetic fields are dragged by galactic circular rotation. We have considered the role of the 
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inclusion of turbulent magnetic diflFusion and the a-effect. If the magnetic turbulent diffusivity is not quenched, the amplitude 
of the warp diminishes in one local diffusivity time-scale ~ H"^ /r)t, which is ?a 0.3 Gyr for rit w 10^® cm^ s~^. Moreover, 
for galaxies moving through the intergalactic medium, the plane-parallel configuration of intergalactic magnetic fields around 
gala^xies, as required in the magnetic model, is by no means expected because of the shear terms. 

Accretion of infalling diffuse material from the intergalactic medium may be inferred from current low-mass star formation 
rates in spiral galaxies. In the model of cosmic infall, the new accreted material reorients the surrounding halo potential, as 
skew angular momentum is absorbed, and antisymmetric integral-sign warps can be excited even if the accreted matter was 
purely dissipative gas (e.g., Jiang & Binney 1999; Bournaud & Combes 2002). The situation is completely different in the case 
of oblique accretion onto the disc of an intergalactic flow caused by the relative motion of the galaxies through the IGM. In 
such scenario, the flow has no net angular momentum at infinite distance but, when it reaches the disc, the impulse transfer 
has an azimuthal dependence and both the m = and the m = 1 modes are excited, generating asymmetric warps. We have 
inferred the amplitude and level of warp asymmetry in the accretion model following a completely different formulation than 
in LC. It turns out that, as reported by LC, an accretion flow velocity of 100 km s~^ and densities of 6 x 10~^ cm~^ can 
account for the amplitude of the Galactic warp. However, the model presents difficulties in explaining the existence of fairly 
antisymmetric warps and, hence, accretion of a plane-parallel flow cannot be the only cause of warps. For the Galactic model, 
the warp should be of type I with the maximum gas displacements on the two sides in the ratio 3 : 2. In addition, it predicts 
that a fraction ~ 3.5% of warped galaxies, at least, should be of type L or U; this fraction may be larger if the clouds have 
a no-null velocity dispersion. There are other predictions that are not borne out by the observations: (1) this model predicts 
that antisymmetric warps should be of type I, whereas those warps falling back towards the midplane of the inner disc should 
present a U-rim shape; (2) if accretion is the only mechanism to excite warps, a strong correlation between amplitude of the 
warp and size of the Hi disc is expected, (3) the rather symmetric vertical distribution of the volume density in the warp 
of our Galaxy indicates that the pressure by the action of the clouds rain on one cap of the H i disc does not produce an 
appreciable vertical distortion; (4) the external pressure required to bend the outer disc would produce H I scaleheights too 
small; (5) the thickness of the Hi disc should have a strong dependence on azimuth. We suggest that the dynamical pressure 
needed to take the gas off the plane is a factor ~ 10 much larger than observations permit. 
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APPENDIX A: ANALYSING THE DYNAMICS OF A DISC UNDER A TORQUE OF TYPE A 

Al Response of the disc and temporal evolution 

We start examining the motion of a disc cloud in Cartesian coordinates {xc, Vc, Zc), originally located in the xy plane (i.e. Zc = 
ic = at t = 0), under the action of an axisymmetric gravitational potential created by the disc and the dark halo, ^{R, z), 
and a vertical nongravitational force /ng {R, <j>) = fi (R) sin (j) as described in case A of §2. Let us consider only displacements 
of the cloud along the z cixis. The equation of motion for the three coordinates are: 

Xc = —fi'xc, (Al) 

j/e = -n^yc, (A2) 

«c = -I'tZc + fi {R) sin (j). (A3) 

Integrating the equations for Xc and j/c, the azimuthal angle of a cloud evolves according to (j) = 4>o + ^t, where 4>o is the 
azimuthal angle of the cloud at t = 0. Hence yc{t; (po) = Rsinijlt + 4>o) and Eq. (A3) becomes 

Zc = -h'tZc + fisin{nt + <j)o)- (A4) 
The solution of this equation after imposing the assumed vertical displacement and speed at t = is 

Zc{t; <j>o) = Zi sm{Qt + (j)o) - Zh sm{utt + (pn), (A5) 
where 



Zh{4>o) = Zl — {cos (j)o) 



2 

2 

1 + arctan 



1/2 



(A7) 



and 

<pH{<t>o) = arctan tan^o^ . (A8) 

As discussed at length in §2, the component zh sm.{i'tt + 4'H) is erased by the winding process and hence, after a few dynamical 
time-scales, only the stationary configuration 

«c(i; ^o) = ^Jl^2 sin(f2* + 0o), (A9) 

is observable. In this particular configuration, a cloud oscillates with frequency Q. in the vertical direction. The vertical 
displacement h{R,(f),t) can be immediately calculated from Zc{t;4>o) by using the identity <po = <p — fit, i.e. h{R,4>,t) = 

Zcit;(l>-nt). 



A2 Balance of torques in the stationary configuration 

In §2 it was shown that 

h{R,<t>.t)= ^J^^^ sin0, (AlO) 

is the only configuration that describes a stationary warped disc in case A. We wish to check that the total torque acting on 
each ring is zero. 

The total (gravitational plus non-gravitational) force per unit mass acting on a parcel of gas at the position (x, y, z) is 
/ = -n^xcx - -I- (/i 8\n4>- i'tz)ez- (All) 

The torque associated with this force is 

T = r X f ^ [{fl^ - ul)yz + yh{R) sin^] 

+ - n^)xz - xfi{R) sin0] By. (A12) 
The torque acting on a ring of radius R is obtained by integrating over the azimuthal angle 

/'27r 

THng(-R)= / Td<j}. (A13) 
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In cylindrical coordinates we know that x = Rcoscj) and y = Rsincj). From Eq. (AlO) it holds that z = {vf — 0^)~^/i sin</). 
Substituting these expressions into Eqs (A12)-(A13) and performing the 0-integration we obtain rring(J?) = 0, implying that 
the total (disc+halo+non-gravitational torque) torque in a ring vanishes, which is fully consistent with the fact that the rings 
do not precess. 



APPENDIX B: EVOLUTION OF THE MAGNETIC FIELD IN A ROTATING DISC WITH VERTICAL 
MOTIONS 

Bl Induction equation and magnetic evolution in a imposed velocity field 

Vertical motions associated with the vertical displacement of the disc and differential rotation can both modify the configu- 
ration of the magnetic lines in the disc. In this appendix we consider the evolution of the magnetic field under field freezing 
conditions in a steady tilted disc. We start with the full equations for dBa/dt, dB^/dt and dBzjdt that result from the 
induction equation (20) by a imposed velocity field of the form 

V = v^{R)e^ + Vz{R,(i>)ez. (Bl) 

In the absence of turbulent motions, the resulting equations are 
BBr 



Dt 



0, (B2) 



where D/Dt = {d/dt + Q,djd4> + Vzdjdz). In the derivation of these equations we have adopted the simplifying assumption 
that dBz/dz « across the thickness of the disc. Since the magnetic force is given by the term BbOBz /dR, we are mainly 
interested in the evolution of Br and Bz. According to Eq. (B2), Br, is constant in an element volume moving with the fluid. 
However, vertical motions could lead to variations of Bz due to the gradients dvz/dR and dvz/d4>. 

In order to estimate the importance of these terms, let us assume that the dynamics of the disc can be described by 
h{R,4>,t) = hi{R)sm{<j} — ujt), where iL>{R) is a known azimuthal frequency which may depend on R. With the help of the 
definition of Vz given in Eq. (12), substituting into Eq. (B4), and exploiting the relations (B2) and (B3), Bz{R,4>,t) can be 
obtained after integration. Define the initial magnetic field B^°^ by 

B^°\R,^) = B{R,<p,0), (B5) 
then Bz is given by the sum of three terms 

Bz (R, 4>,t) = Bf^ {R, 4>-nt) + B<"' {R, 4>,t) + Bi"^ {R, 4>, t) + B<'=' {R, 4>, t) , (B6) 
where 

B'f ^ {R, <!>, t) = Bf> {R, <l>-^t)^ [sin (0 - ut) - sin - Ot)] , (B7) 

B*'"' {R, 4>, t) = Bf ( J^, - m) ^ [cos (0 - ujt) - cos (0 - m)] , (B8) 
and 

B'f^ {R, <p, t) = Bf {R, (t> - Qt)hi ^(^Iz^t cos (0 - ujt) . (B9) 

Although this solution is mathematically exact, the growth of magnetic gradients will perturb the fluid motions v and, 

therefore, this solution is only valid at times short enough that one can ignore the back reaction. 

The first special case occurs when u!{R) = n{R). Under this condition M"' = M*' = Bi"^ = and therefore Bz{R,(t>,t) = 
Bi°^ {R, 4> — Q,t). The question arises as to how dBz/dR varies. To do so, one has to remind that 

dB(°\RA-^t) _ dB(°\R,4>) dB(°\RA) dQ^ 

dR ~ OR d4> dR ^ ' 

Therefore, if initially Bz depends on the azimuthal variable <j), then dBz/dR grows linearly in time due to differential rotation. 
dBz/dR will be time-independent only \i u) = Q. and either Bz is axisymmetric or the disc rotates in a solid-body manner 
(or both). Summing up, in a differentially rotating disc with u! = Q. and Bz being a function only of R, the forcing term 
BndBz/dR is constant in an element volume moving with the fluid. In the following we discuss the case u> ^Q,. 
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In the magnetic scenario for warps formation, we know from Eq. (19) that both _B]j' and ' depend on the azimuthal 
angle <p. Thus, the radial derivatives of B^°'\ B^^ and 5^°' will contain a term oc tdQ/dR due to the second term of the RHS 
of Eq. (BIO). Hence, the magnetic gradient dBz/dR will grow in time and, therefore, no steady-state configuration exists in 
a differential rotation disc. 

In addition, from Eqs. (B7)-(B9) it follows that, in the presence of vertical motions^ {w ^ Q), Bz and hence dB^/dR 
depend on (p, even if Bz were initially independent of (f>. Therefore, the magnetic term BadBz/dR, which was supposed initially 
of the form sin((;6 + C), with Q a constant phase, as required to drive integral-sign warps (see §3), will develop undesired modes 
with m> 1, for example, oc sin((/)-|-^) sin(0-|-f) and so on. These modes of higher wavenumber would spoil the original m = 1 
shape of the warp. In the next Subsection we will estimate the evolutionary time-scales for the growth of magnetic gradients. 

Finally, we must remark that Muj ^ Q, then Bz{R, 4>, t) will depend on Hence, the term {pR)~^B^dBz/d4> that appears 
in Eq. (18) increases with time. 



B2 Evaluation of the magnetic evolutionary time-scales in a galactic disc 

As discussed in the previous Subsection, in order to preserve the ability of the magnetic torque to maintain a m = 1 distortion 
during a time ta, the magnetic gradients should change only slightly in that time. We introduce the characteristic growth 
time-scales of the magnetic stress component BrOBz/OR, r'"' with n = a,b,c, as the times that are required for dBi"^/dR 
to change by of order of dBf^ /dR itself, according to the following relation 



dR 



''"^ (BU) 



We are now in a position to estimate these characteristic time-scales by performing the derivatives with respect to R in 
Eqs. (B7)-(B9). 

Given x(^) & certain variable that depends on 7?, it is convenient for notation to define Ax s \x[R2) — x(^i)l) where 
R\ and R2 are the minimum and maximum radii of the warp. For example, A/ii = |fti(7?2) — h-i(Ri)\ = hi{R2) is the height 
of the warp at i?2. With this notation, the tip warp angle denoted by dm is given by tan dm = hi_{R2) / {R2 — Ri) = Ahi/AR. 

The characteristic growth time-scales for the different contributions are 

j(0) 



(a) 

' = max 



ABl^AR_i_ n , 

^ 5(0) A/ii An' 2\n-u\^' ^ ' 

where the bar denotes a characteristic value in the warp, say x = x([-R2 -|- Ri\/2). Analogously, 
(b) AbI[AR 1 TV , 

and 

r<"' =min{Ti,T2,r3}, (B14) 
with 

ABi°^ AR 1 

(abPar 1 

and 



/ ABi°^ AR 



1/2 



In order to give guide numbers of these time-scales, let us explore the limit case in which the amplitude of the warp is 
sufficiently high, h ^ (ATTUt p)^^ Br^iOBz/OR, that the magnetic forces do not affect the vertical motions. In such case we 
know from §2 that lo = Q — Vf To estimate representative values for Afl, Auj and Auj' , we employ the same model than 
Binney (1992). Assuming a fiat rotation curve and an exponential disc with scale length Ra, to{R) can be approximated by 

u>{R)c^nf^yQ{R), (B18) 



^ Remind that Vz = {d/dt + Q,d/d(j>)h. Consequently, Vz ^ implies oj ^ CI. 
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beyond ARd- For Ri = 4Rd and R2 — 6Rd, it follows that Aw ~ 0.4Af2 and Alo' ~ Af2. Regarding AS", the following upper 
bounds can be readily derived using the framework of Battaner & Jimenez- Vicente (1998) for magnetically-generated warps: 
ABi°VAB^°' « AB^V-B^j"' « ABi°^/Bf^ < 4. The resulting upper values for a galaxy like the Milky Way with Rd w 3.5 

kpc and Vc ~ 220 km s~^ are 

r'") and t'*'' < Gyr, (B19) 

tan^u, 



r 



< Gyr. (B20) 



(0 ^ 0-6 
Vtan( 

Inserting a typical value for tan^uj of ~ 0.2, a conservative upper value for the evolutionary time-scale is < 1.3 Gyr. 

In order for BudB^ / dR to preserve its initial wavonumbcr m — 1, uj(R) must differ only slightly from fl{R) at every 
radii, i.e. Vz should be very small or zero. The limit Vz —^ corresponds to what we call the "particular solution". When Vz 
is small, the induction equation simplifies to Eqs. (21)-(23). 



